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Abstract 

It is now a firmly established fact that all family lepton numbers are violated 
in Nature. In this paper we discuss the implications of this observation for 
future searches for rare tau decays in the supersymmetric see-saw model. Using 
the two loop renormalization group evolution of the soft terms and the Yukawa 
couplings we show that there exists a lower bound on the rate of the rare process 
/i ^ 67 of the form BR{^ 67) ^ C x BR{t fi'y)BR{T 67), where C is 
a constant that depends on supersymmetric parameters. Our only assumption 
is the absence of cancellations among the high-energy see-saw parameters. We 
also discuss the implications of this bound for future searches for rare tau decays. 
In particular, for large regions of the mSUGRA parameter space, we show that 
present i?-factories could discover either r — > /U7 or r — > ej, but not both. 
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1 Introduction 



The renormalizable part of the Standard Model Lagrangian is invariant under four 
global f/(l) symmetries, namely baryon number, B, and the three family lepton num- 
bers, Le, Ln and Lr. This invariance has been for many years considered acciden- 
tal and expected to be broken by additional terms in the Lagrangian, possibly non- 
renormalizable. Whereas experiments searching for proton decay have not provided yet 
any evidence for baryon number violation, it is nowadays a firmly established experi- 
mental fact that the three family lepton numbers are violated in the neutrino sector. 
Namely, the disappearance of electron neutrinos on their way from the Sun indicated by 
the Homestake chlorine detector [1], SAGE [2], GALLEX/GNO [3,4], Kamiokande [5], 
Super-Kamiokande [6] and Borexino [7], and unequivocally confirmed by SNO [8], 
proves the violation of L^,. This is further supported by the disappearance of electron 
antineutrinos observed by the reactor experiment KamLAND [9]. On the other hand, 
the atmospheric neutrino anomaly discovered by Kamiokande and 1MB [10], and ex- 
plained by Super-Kamiokande [11], Soudan2 [12] and MACRO [13] as an oscillation 
of muon neutrinos into a different neutrino species, proves the violation of L^. The 
disappearance of muon neutrinos reported by the long baseline accelerator experiments 
K2K [14] and MINOS [15] supports this conclusion. Finally, the observation of tau neu- 
trino appearance in the atmospheric neutrino flux by Super-Kamiokande [16] indicates 
the violation of L^. 

The most economical way to accommodate the family lepton number violation in the 
Standard Model is by adding to the leptonic Lagrangian a dimension five operator [17] 

Acp = -e%^Y,,^L,H* - -f{UH){L,H) + h.c. , (1) 

where Li and are the left-handed and right-handed leptons respectively, i^j = 1, 2, 3, 
H is the Higgs doublet and A a mass parameter. In fact, this operator not only violates 
the three family lepton numbers but also the total lepton number. After the electroweak 
symmetry breaking, a Majorana mass term is generated for the left-handed neutrinos, 
A4ij = ^{H'^y, being the smallness of neutrino masses attributed to small values of 
the couplings aij and/or to a large value of A. 

In this minimal framework, lepton flavour violation is generated by the same operator 
that generates neutrino masses. Therefore, the rate of any lepton flavour violating pro- 
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present bound 


projected bound 


BR{fx 67) 
BR{t 67) 
BR{t fi-f) 


1.2 X 10"^^ [19] 
1.1 X 10"^ [21] 
6.8 X 10-s [23] 


10"^^ [20] 
10-9 [22] 
10-9 [22] 



Table 1: Present and projected bounds on the rare lepton decays. 

cess will be proportional to the neutrino masses. Indeed, a detailed calculation shows 
that [18]: 

2 

BR{lj ^ ku^Ui) , (2) 

which gives BR{t fi-f) ~ IQ-^^, BR{fi 67) ~ 10-^^ BR{t 67) ~ 10-^^ far 
below the sensitivity of present and projected experiments (see Table 1)^. 

Nevertheless, the Lagrangian Eq. (1) describes just an effective theory and new 
degrees of freedom are expected to arise above the scale A. The interactions of the 
new degrees of freedom with the lepton fields are likely to contain additional sources 
of flavour violation that can enhance the branching ratios of the rare decays by many 
orders of magnitude, bringing them to the reach of future experiments. For this reason, 
rare lepton decays are considered very powerful probes for physics beyond the Standard 
Model. 

The supersymmetric (SUSY) see-saw mechanism is probably the best motivated 
high energy theory to generate small neutrino masses [25] . In this framework the 
particle content of the Minimal Supersymmetric Standard Model (MSSM) is extended 
with three right-handed neutrino superfields, ujn, i = 1,2,3, singlets under the Stan- 
dard Model gauge group. Imposing i?-parity conservation, the leptonic superpotential 
reads: 

W^iep = e%^Y,,^L,Ha + i/^,Y,,,L,i/„ - ^i^^,M,,z/^ . , (3) 

where and are the hypercharge +1/2 and —1/2 Higgs doublets, respectively, Yg 

and Yi, are the matrices of charged lepton and neutrino Yukawa couplings, respectively, 

and M is a 3 X 3 Majorana mass matrix. It is natural to assume that the overall scale 

^In the Table we restrict ourselves to bounds that were published by the time of writing this paper. 
We note however that the Belle Collaboration has reported the bounds BR{t nj) < 4.5 x 10"^ 
and BR{t ej) < 1.2 x 10"^ in the yet unpublished preprint [24]. 
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of M, which we will denote by Mmaj, is much larger than the electroweak scale or any 
soft mass. If this is the case, at energies below Mmaj the theory is well described by 
the following effective superpotential: 

W^iS = 4.Y,,,.L,if, + i (YjM-^Y,)^^^. {UH^){L,H^) , (4) 

that generates the fermionic Lagrangian Eq. (1). In the phenomenological studies it is 
convenient to work in the leptonic basis where the charged lepton Yukawa coupling is 
diagonal, Yg = diag(2/e, Uiii Ut)- Then, in this basis, the neutrino mass matrix is given 

by 

M = (YjM-^Y,) {H'^f , (5) 

whose eigenvalues are naturally very small due to the suppression by the large right- 
handed neutrino mass scale. 

From the point of view of generating small neutrino masses, the non-supersymmetric 
version of the see-saw mechanism is equally natural. Nonetheless, in the non-super- 
symmetric see-saw model the presence of very heavy new degrees of freedom interacting 
with the Higgs doublet introduces a serious naturalness problem. The Higgs mass ac- 
quires quadratically-divergent radiative corrections which would drive the Higgs mass 
to values of the order of the Majorana mass scale [26]. However, in the supersymmetric 
version this divergence is automatically canceled by the presence of right-handed sneu- 
trinos with a mass essentially identical to the mass of their corresponding right-handed 
neutrinos. Therefore, the supersymmetric see-saw mechanism can accommodate simul- 
taneously tiny neutrino masses and a relatively light electroweak scale without serious 
fine-tunings. 

In addition to the flavour violation stemming from the right-handed neutrino Yukawa 
couplings, the supersymmetric see-saw model contains additional sources of lepton 
flavour violation in the soft SUSY breaking Lagrangian [27]: 

(^AeijFfiiHdLj + A^ijiy*j^^HuLj + h.c. j + etc . (6) 

where Lj, ejn and are the supersymmetric partners of the left-handed lepton dou- 
blets, right-handed charged leptons and right-handed neutrinos, respectively, m|, 
and are their corresponding soft mass matrices squared, and Ag and Aj, are the 
charged lepton and neutrino soft trilinear terms. 
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The flavour violation in the slepton sector contributes through one loop diagrams to 
different flavour violating processes such as rare muon and tau decays, — > e^/i^ or 
yU — e conversion in nuclei. The strong bounds on these processes restrict very severely 
the structure of the soft mass matrices at low energies [28] , suggesting an approximately 
flavour universal structure: (m|)jj ~ rrij^Sij, ~ rnl6ij, (Ae)ij — AeYeij. The 

most plausible explanation for this structure is to assume that supersymmetry is broken 
in a hidden sector and the breaking is transmitted to the visible sector by a flavour 
blind mediation mechanism. If this is the case, the soft terms would be strictly flavour 
universal at some high energy scale: 



Interestingly, if this high energy scale is larger than the right-handed neutrino 
masses, the flavour violation in the neutrino Yukawa couplings will propagate through 
radiative effects to the soft terms [29] . As a consequence, even under the most conserva- 
tive assumption for the soft terms from the point of view of lepton flavour violation, in 
the super symmetric see-saw model some amount of flavour violation is always expected 
at low energies. 

The discovery of small neutrino masses as a hint of the see-saw mechanism and 
the continuous improvement in sensitivity of the experiments searching for rare lepton 
decays have stimulated in recent years the interest in the radiative generation of lepton 
flavour violation. Different groups have estimated the predictions for various lepton 
flavour violating processes in some speciflc high-energy frameworks based on Grand 
Unification [30], flavour models [31], texture zeros [32] or string theory [33], but also 
pursuing a more phenomenological approach, where the predictions are somehow con- 
nected to low energy observables [34-38]. A general expectation of all these analyses 
is that the observation of rare lepton decays could be possible at the next round of 
experiments. Nevertheless, it should be stressed that there are many optimistic as- 
sumptions underlying this conclusion, and the non-observation of rare decays in future 
experiments would by no means exclude the supersymmetric see-saw model. 

In this paper we will adopt a completely phenomenological approach and we will de- 
rive the relation among the branching ratios BR{fi 67) ^ C x BR{t —* fi^)BR{T 
67) that holds independently of the high-energy see-saw parameters, with the only 








(7) 
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assumption of the absence of cancellations. In this relation, however, some model de- 
pendence will remain associated to our ignorance on the supersymmetric parameters, 
which is encoded in the constant C. To derive our main result we will present in Sec- 
tion 2 the scenario that saturates the bound and yields the minimal rate for — > 67. 
In Section 3 we will compute the off-diagonal terms of the slepton mass matrices in 
this extremal scenario taking into account the two loop renormalization group equa- 
tions and we will derive relations among them. Using these relations, we will derive in 
Section 4 the lower bound on BR{fi — > 67) in terms of the branching ratios of the rare 
tau decays. We will also discuss the implications of this bound for present and future 
searches for rare tau decays for different supersymmetric benchmark points. Finally, in 
Section 5 we will present our conclusions. We will also present an Appendix containing 
the complete two loop renormalization group equations for the soft SUSY breaking 
terms, including the effects of the neutrino Yukawa couplings, which to the best of our 
knowledge have not been given explicitely in the literature. 

2 Scenarios with minimal rates for /i ^ 67 

In the MSSM extended with right-handed neutrinos there are sources of lepton flavour 

violation both in the soft SUSY breaking Lagrangian and in the SUSY conserving 

Lagrangian. The latter can be completely encoded in the neutrino Yukawa couplings 

by means of a basis transformation. Neutrino oscillations require the presence of lepton 

flavour violation in the neutrino Yukawa couplings, but not in the soft SUSY breaking 

Lagrangian. Therefore, the minimal rate for /i — 67 will clearly occur in scenarios 

where there is no lepton flavour violation in the soft SUSY breaking Lagrangian. There 

are in fact some well motivated supersymmetric scenarios where the soft breaking 

parameters are flavour universal at some high energy scale. However, if this scale 

is larger than the mass of the right-handed neutrinos^, the flavour violation in the 

supersymmetric part of the Lagrangian will propagate to the soft SUSY breaking terms 

through quantum effects. 

For generic neutrino Yukawa couplings, the off-diagonal elements of the soft SUSY 

^This is the case for example in scenarios with minimal supergravity, dilaton-dominatcd SUSY 
breaking, or gauge mediation where the masses of the messenger particles are larger than the masses 
of the right-handed neutrinos. 
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breaking terms read at low energies, in the leading-log approximation, 

K),^. - -^K + m^ + <„ + |A.P)(YtY.).,log(^) , 

where ^ 7^ j and A/x is some high energy scale that we identify with the Grand Unifi- 
cation scale, Mx = 2 x 10^^ GeV. The size of the off-diagonal elements depends very 
strongly on the flavour structure of the neutrino Yukawa couplings, which in the see- 
saw model is not directly connected to the flavour structure of the low energy neutrino 
mass matrix. In fact, there is an infinite set of neutrino Yukawa couplings compati- 
ble with a given set of low energy data [36]. Among all those Yukawa couplings, the 
minimal rate for /i — 67 will clearly occur in the scenario where 

(YlY,)^^ (M^^j) = yl^yu + 2/2*22/21 + 2/322/31 = . (9) 

Assuming that there are no cancellations among the different terms, this condition 
is satisfied in the following eight situations: 



(a) 


yn 


= 0, 2/21 


= 0^ 


-2/31 


= 0, 


(e) 


2/12 = 


0, 


2/21 = 


0, 


2/31 = 


0, 


(b) 


yii 


= 0, 2/21 


= 0, 


2/32 


= 0, 


if) 


2/12 = 


0, 


2/21 = 


0, 


2/32 = 


0, 


(c) 


yii 


= 0, 2/22 


= 0, 


2/31 


= 0, 


(g) 


2/12 = 


0, 


2/22 = 


0, 


2/31 = 


0, 


(d) 


yii 


= 0, 2/22 


= 0: 


■ 2/32 


= 0, 


(h) 


2/12 = 


0, 


,2/22 = 


0, 


,2/32 = 


0, 



(10) 

where the Yukawa matrix elements are evaluated at M^aj- The cases (a) and {h) pre- 
serve electron and muon lepton numbers respectively. Since in this basis the neutrino 
Yukawa coupling is the only source of lepton flavour violation, the complete Lagrangian 
will preserve at least one family lepton number. This invariance is inherited by the ef- 
fective theory, in conflict with the flavour conversions observed in neutrino oscillations. 
On the other hand, the remaining six possibilities (6) — (g) violate all family lepton 
numbers and lead to a neutrino mass matrix that schematically reads 



~ X X (11^ 
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and which is allowed by present experiments (this mass matrix leads to the prediction 
sin ^13 ~ 2 V Am^"°' si'^^^soi — 0.08 for the case with normal neutrino mass hierarchy 



and sin — ^^2°°' sin26'soi — 0.04 for the case with inverted neutrino mass hier- 
archy). Therefore, from Eq. (10) it is straightforward to check that the experimentally 
allowed neutrino Yukawa textures which lead to a vanishing rate for /i — ^ 67 necessarily 
lead to (Yi[Y^)3i(Mmaj) 7^ and (Yi!Y^)32(Mmaj) 7^ 0, unless unnatural cancellations 
in Eq. (9) are taking place. 

The presence of exact zeros in the Yukawa matrices as in Eq. (10) can be justified by 
symmetries. However, symmetries are not expected to hold at the decoupling scale but 
at the Grand Unification scale, Mx. Again, and barring cancellations, the minimal rate 
for yU — > 67 will occur when {YI'Y y)2i{Mx) = 0, that as was argued above necessarily 
implies (YtY.)3i(Mx) ^ and (YtY.)32(Mx) 7^ 0. 

The key point of this paper is that even when the soft terms are flavour universal 
at Mx and {Y^Y y)2i{Mx) = 0, the flavour violation in the r — yU and the r — e sectors, 
will propagate through radiative corrections to the /i — e sector, inducing small though 
non- vanishing values for (m|)2i and (Ag) 12.21 at low energies. We will show in this 
paper that in this very special case the rate of /x — 67 is related to the rates of the 
rare tau decays through 

BR{^i ^ 67) ~ C X BR{t fi-f)BR{T 67) , (12) 

where the proportionality constant, C, will be determined in the next sections. 

This value is attained by the experimentally allowed neutrino Yukawa textures 
satisfying (Y;'!Y,^)2i(Mx) = 0, which is the scenario producing the minimal rate for /i 
67. Therefore, and barring cancellations, for any other Yukawa coupling compatible 
with the low energy neutrino parameters the following bound will hold: 



BR{fi ^ 67) > C X BR{t fi-f)BR{T ^ 67) 



(13) 



This inequality is the main result of this paper. As we will see later, this bound 
has important implications for the searches for rare tau decays in present and future 
experiments. 

To finish this section we would like to comment on the possibility sometimes dis- 
cussed in the literature of establishing model independent correlations between the 
branching ratios of the rare decays of the form BR{t ^'-f)/BR{fi 67) ^ constant. 
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We would like to clarify this point stressing that in a general see-saw model with three 
right-handed neutrinos such correlations cannot be established on model independent 
grounds. 

It can be proved that there is a one to one mapping between the high energy see-saw 
parameters, and M, and the matrices A4 and Yj^Y,^ [39]. The former is the neutrino 
mass matrix and the latter, the matrix that participates in the radiative corrections 
to the soft SUSY breaking terms, which is in principle measurable at low energies 
once the boundary conditions for the soft terms have been specified. In particular, 
the off-diagonal elements of YjiYj, are directly related to the branching ratios of the 
rare decays. Being this mapping bijective, one could use as parameters of the see- 
saw model either the familiar set {Y^, M} (the top-down parametrization) or the less 
familiar one {A4,YIY^} (the bottom-up parametrization). Clearly, in the bottom-up 
parametrization the branching ratios for the rare decays are inputs, and as such are 
completely uncorrelated. For the same reason, it is not possible to establish on model 
independent grounds correlations between the branching ratios of the rare decays and 
neutrino parameters, such as sin ^13, the neutrino mass spectrum or the CP violating 
phases. This result also holds for the most restricted case with just two right-handed 
neutrinos [40]. 

A more explicit proof of this result can be derived as follows using the familiar top- 
down parametrization. Let us first write the effective neutrino mass matrix in terms 
of the high-energy parameters in the basis where the charged lepton Yukawa coupling 
and the right-handed mass matrix are simultaneously diagonal: 

One could conceive a scenario where -C ^^j^, ^^j^ for all i,j = 1, 2, 3. This could 
occur for instance when the mass of the heaviest right handed neutrino mass is much 
larger than the masses of the other two right-handed neutrinos. If this is the case, the 
third row of the neutrino Yukawa matrix is completely unconstrained from neutrino 
observations. In addition, it could occur that y^i are much larger than yn, y2i, implying 
that (YlY,^)ij ~ y^iysj- Therefore, the branching ratios for the rare decays would be 
essentially determined by the parameters y^i, that are completely unconstrained from 
low energy data. As a result, the ratio BR{t fi'-f)/BR{fi — > 67) can range from 
zero (when 2/33 = 0, 1/31,2/32 7^ 0) to infinity (when y^i = 0, 1/32,2/33 0), being both 
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situations compatible with neutrino observations. Additionally, since the parameters 
that determine the branching ratios, ?/3j, are unconstrained from low energy data, no 
model independent correlation can be established between rare decays and neutrino 
parameters. 

In order to establish correlations between two branching ratios, or correlations 
between one branching ratio and low energy neutrino parameters, it is necessary to 
make assumptions about the high-energy see-saw parameters, for instance motivated 
by Grand Unified theories or by flavour models. In this paper we make what we consider 
the most minimal assumption about the high-energy see-saw parameters, namely that 
in constructing the derived quantities A4 and YIY^ no cancellations occur. Then, 
the no-go theorem presented in [39,40] can be circumvented, opening the possibility 
of deriving relations among the branching ratios. Although the relation that results, 
BR{fi — > 67) ^ C X BR{t fi'y)BR{T — > 67), is not model independent in the 
strict sense, it is valid in a very large class of models, characterized by the absence of 
unnatural cancellations. In addition, this bound will proof to have remarkably strong 
implications given the generality of our assumptions, as we will see in Section 4. 

3 Radiative generation of lepton flavour violation 
in the soft mass matrices 

The Renormalization Group Equations (RGEs) for the relevant soft terms read 

'""^ - \/32+7T^/3S+.... (15) 



L 



dt levr^^^^i (167r 

dt ~ 167r2'^-i+(167r2)2'^-i + -' 

dt ~ 167r2'^^^+ (167r2)2^A« + -' > 

where (3'^^'^ and indicate respectively the one and two loop /3-functions. The com- 
plete expression for the /3-functions of the soft terms in the supersymmetric see-saw 
model can be found in the Appendix. 

The solution to the RGEs can be well approximated by the trapezium rule, which 
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keeping just the one and two loop contributions reads: 

1 1 



maj; 



mi(Mx) - 



1 



maj; 



Ae(M 



maj; 



(167r2; 


22 


m2(Mx; 


- 


1 


1 


(167r2; 


22 


Ae(A'/x) 




1 


1 



167r2 2 I '"i 

" 5(2) 



maj; 



1 1 



maj; 



/5;^UMx) + /5«(M^aj; 



log 



log 

Mx 



M- 



X 



M 



maj 



M, 



maj 



167r2 2 



log 



log 

Mx 



Mv 



M 



maj 



M 



maj 



(19) 



1 1 



167r2 2 



majy 



(167r2)2 2 



/3g(Mx)+/3i^i(M^aj) log 



(2), 



log 

M 



M 



X 



M, 



maj 



X 



M, 



maj 



(20) 



For the Yukawa couplings leading to a minimal rate for /i —* e-f, Eq. (10), the 31 
and 32 entries in the left-handed slepton mass matrix are generated at order 0{Y^). 
Therefore, noting that (3^^l (M^aj) = P^^l (Mx) + OiY^) and that the two loop /?- 



functions are 0{Y^), it follows from Eq. (18) that 



mi(M, 



maj j31,32 



167r^ 



(/52)31,32(Mx)log 



M 



X 



M, 



maj 



(21) 



Using Eq. (75) in the Appendix for the one loop /9-function we find 



mi(M, 



maj j 31,32 



57r 



:{ml + ml + ml^ + \A,\^) (Y^Y,) 31,32 log| 



M. 



M 



maj. 



+ 0{Y,') , (22) 



where all the soft terms in the right-hand side of this equation and (Yi^Yi,)3i,32 are 
understood to be evaluated at the scale Mx- 

On the other hand, since (Yi^Yi,)2i(Mx) = the previous equation indicates that 
m2 (Mmaj)2i vauishes at order 0(Y^) and is only generated at higher order in per- 
turbation theory. Keeping only terms of (9(yj'), it follows that {l3^2)2i{M„ 
(/^S hi{Mx) + 0{Y^) and (/^^J^ )2i(Mx) = 0. Therefore, 



'majj 



>i)2i(M, 



maj; 



^o(^m?)2l(^maj) 



167r2 2 ''^'"t 



log' 



M 



X 



M 



maj 



'167r 



^(/52.)2i(Mx)lo^ 



M. 



M 



+OiY^ 



maj/ 

(23) 



The one loop /3-function at the Majorana mass scale is proportional to (YjiY,^)2i(Mi 



maj / ) 



that following Eq. (56) in the Appendix is generated radiatively by the effect of 
(Y^Y,)32 and (Y1Y,)3i. The result is: 



(YtY,)2i(M^,j) 



i(YtY.);2(YtY.)3ilogf 



57r 



M, 



maj 



+ 0{Y^) . 



(24) 
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Using Eqs. (23,24) and Eqs. (75,76) in the Appendix, it is straightforward to com- 
pute (m|)2i(Mmaj)- The result is: 

1 



(mi)2i(M, 



maj> 



X 



121og| 



167r2)2 

2 



M- 



{mi + mjj^ + mt + 2\A, 
Mx 



+ 8 logi 



YiYi,)32l 



YtY,U + OiY,') . (25) 



'maj/ V"maj/ 

Recall from the previous section that when (YIY^)2i{Mx) = 0, the observed pattern 
of neutrino mixing angles requires {YlY^)'ii{Mx) 7^ and {YIy^)32{Mx) 7^ 0. Then, 
(m|)2i(Mmaj) will always be generated at least at order 0(Yjf): it will be generated at 
order 0(Yj^) when (YIY^)2i{Mx) is different from zero and at order 0{Y^) when it is 
equal to zero. Therefore, and barring cancellations, the observed violation of all family 
lepton numbers in neutrino oscillations necessarily implies in the supersymmetric see- 
saw model a contribution to the process fi ^ c-f from the off-diagonal entries of the 
soft terms. 

Following the same steps , one can calculate the off-diagonal elements of the trilinear 
soft terms. For the elements generated at 0{Y^) we obtain 



(A, 



1 



6^31,32 



(A, 



ej 13,23 



167r2 
1 

167r2 



logI 



logI 



X 



A^maj 
Mx 



M, 



maj 



(2A, 
(2A, 



Ae)y.(YtY,)3l,32 
Ae)ye,^l{YlY,y)'^l 



32 



0{Y, 



(26) 
(27) 



and for the elements generated at 0{Y^ 

2 



7 / Mv 



(A. 



167r2 

xy,,(YtY,)*2(YtY.)3i 
7 



maj 



e 12 



167r2 



(14A^ + -Ae)log 



0{Y,') , 

Mx 



'maj 



+ + 2A,) logI 



[SA, + 2Ae) logI 



X 



M, 



maj 



X 



M 



maj 



xy,(YtY.)*i(YtY.)32 + 0(i;^ 



(2^ 



(29) 



Finally, there are additional sources of lepton flavour violation stemming from 
the radiatively generated off-diagonal entries in the charged lepton Yukawa coupling. 
Analogously to the previous discussion, (Ye)(Mmaj)i2 and (Ye)(Mmaj)2i get values of 
OiYj"), while the remaining off-diagonal terms get values of 0{Y^). Rotating the 
leptonic fields to bring them to the basis where the charged lepton Yukawa cou- 
pling is diagonal will modify the values of the soft terms calculated above. Defining 
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Ye = VnYf^'Wl where Yf^s is a diagonal real matrix and Vl,r are unitary matrices, 
the basis transformation L VlL, cr — >■ Vnen yields 



m, 



(30) 



The explicit expression for Vl is 



Vr 



1 

_T/* _T/* 




(31) 



where 
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5log(^ 
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167r2 



5, fM^ 
- log 

2 ^VM, 
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21og| 



X 



^^maj 
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L, 13.23 



IGvr^ 
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maj 



(YtY,)*,(YtY,)32 + 0(i;^), 

YtY,)*,(YtY,)3i + 0{Y^') , 

(32) 



Notice that Vl,2i 7^ ""■^l,i2) ^.s required by unitarity. On the other hand, the expression 
for Vr is 



V, 



R 



1 Vr,12 yR,13 

V/J,21 1 Vr^23 

— V* —V* 1 

^R,13 ^R,23 ^ 



(33) 



where 

VrA2 



Ve 



VrJ V^majy 



logi 



X 



maj 



[YtY,)*,(YtY,)32 



ii,21 



V, 



1 2/o 



ye 



logI 



X 



maj 



logI 



M- 



X 



ii,13 



log' 



X 



iJ,23 



1 I/m, /^x 



(^^^1^)31 



(Y,iY,y)32 



'maj 

r4\ 



(YtY,)*2(YtY.)3i 



(34) 



With these expressions for Vl and Vr it is straightforward to compute, using 
Eq. (30), the off-diagonal elements of the leptonic soft mass terms in the basis where 
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the charged lepton Yukawa couphng is diagonal. In this "mass basis" the off-diagonal 
elements of the right-handed slepton mass matrix read: 



m 



.2\mb 
eJ31,32 



(35) 



that no longer vanish, but still give negligible contributions to the rare decays compared 
to the other sources of flavour violation. On the other hand, for the left-handed slepton 
mass matrix they approximately read 



maj; 



1 



167r^ 



{ml + ml+ml + 2\A,\^)\og\ 



M, 



maj 



{ml + +ml + ^\A,\'') log(^ 



maj 



(YtY.)*2(YtY,)3i 



(36) 



Note that the diagonalization of the charged-lepton Yukawa coupling does not alter 
m|(Minaj)3i,32 but ouly (m|)(Minaj)2i- Finally, the off-diagonal elements of the soft 
trilinear term read: 



^mb 
^ej3l,32 



J 13,23 



1 , /Mx 
log 



(A, 



mb 
eJl2 



1 

167r2 



167r2 



logi 



Mx 



M 



maj 



2A,y,(YtY,)3i,32 



2yljy/le,^(Y;J^Y,y)3-|^ 32 , 



81og| 



81og| 



X 



M ■ 

-"-'maj 

Mx 



maj 



Slog 



Slog 



M, 



maj 



M, 



maj 



A,y^(YtY,)*2(YtY,)3i , 
A,y,{YlY,)l,{YlY,),2 • (37) 



In this case, all the off-diagonal elements get modified at the lowest order by the basis 
t r ansf ormation . 

It is apparent from Eqs. (36,37) that in the scenario with flavour universality of 
the soft terms at Mx and (YjiYj,)2i(Mx) = there exists a very precise correlation 
between the 21 and the 31 and 32 entries of the soft terms. As was argued in the 
previous section, this scenario yields the minimal amount of flavour violation in the 
fi — e sector. Therefore, any other neutrino Yukawa coupling compatible with the 
experimental data will induce larger (m|)^^ and (Ae)5^'^2i entries at low energies. For 
example, under the assumption of complete universality of the soft mass terms at 
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Mx, m|(Mx) = m|^^(Mx) = ml{Mx) = rnl, the following lower bounds hold at low 
energies for the 21 entries of the soft terms: 

2 



m 



2 \mb 
L)21 



> 



2{3ml + 2|A.n log(-g^) + 2{3ml + ||A.n log(^ 



(A \"^^( \ 



(3m2+|A,|2)log| 
2 



Mx 

Mm a 



> 



logl 



Mx. 



(3^ 



These bounds will eventually translate into a bound on BR{fi 67) involving the 
branching ratios of the rare tau decays. 

4 Lower bound on fi ^ from rare tau decays 

After computing the radiatively generated off-diagonal elements of the soft SUSY 
breaking terms, it is straightforward to calculate the branching ratios for the rare 
lepton decays. In order to understand qualitatively the results we will use in this sec- 
tion approximate formulas for the branching ratios, although in our numerical analysis 
we used the general expressions existing in the literature [34] and we solved numerically 
the renormalization group equations including the two-loop /3-functions. 

A very useful tool to treat analytically the complicated exact expressions for the 
branching ratios is the mass insertion approximation, where the small off-diagonal 
elements of the soft terms are treated as insertions in the sfermion propagators in the 
loops [28,35]. This rationale can also be applied to the gaugino-higgsino sector, yielding 
at the end of the day relatively compact expressions for the branching ratios. 

It was argued in [35] that for the rare tau decays, the dominant contributions 
correspond to the mass insertion diagrams enhanced by tan /? factors. Namely, 

l(mi)3iP 



BR{t 
BR{t - 



tan f3BR{T ei^rZ/g) 



l(mi)32p 



/^2 



rrvi 



tan^ [3BR{t — > fii'r^^) 



(39) 



where BR{t fiUr^n) — 0.17, BR{t eVrVe) — 0.18 and ms is a mass scale of 
the order of typical SUSY masses. In the case of the Constrained MSSM it is best 
approximated by m| ~ {).'bm'^My^{rn^ + 0.6M^y2)^ [38], where mo is the universal 
scalar mass and M1/2 is the universal gaugino mass at Mx. 
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On the other hand, being |(m|)2i| and |(Ae)2i| generated at order 0{Y^), we keep 
for consistency the contributions to BR{fi — > 67) induced not only by a single mass 
insertion (where (m|)2i or (Ae)2i is inserted) but also by a double mass insertion 
(where {m\)s2 or (Ae)32, and (m|)3i or (Ae)i3 are inserted). The result reads: 



a' 



BR{fi 67) _ ^2 



[mlhi (mi)*2(mi)3i 



tan^ p , (40) 



where ms' is another mass scale of the order of m5, although in general different 
(note that the single and the double mass insertion contributions have opposite signs). 
Inserting the bound Eq. (38) into Eq. (40), we obtain the following bound: 

M(M-e,)>^MkQ(EiMw^, (41) 

where rris" is another mass scale, again of the same order of ms, ttls'. Using the 
expressions for the rare tau decays in terms of |(m|^)3ip, |(m|)32p, Eq. (39), this 
bound can be casted as: 

BRU e7) > al. BR[r^^^,) BRjr ^ e,) 

'^^ ~ a3 tan2 {3 mf„ BR{t /iz/,P^) BR{t ^ evr^e) ' ^ ' 

This equation is a more explicit expression of Eq. (13) and constitutes the main result 
of this paper^. 

The numerical values of and rris" depend crucially on the supersymmetric sce- 
nario. Before presenting the exact numerical results for some common supersymmetric 
scenarios, let us obtain first a rough estimate of the bound Eq. (42). To this end, we 
will make the approximation ms" ~ ms- Then, the previous bound reads: 

Rfi-r ^ > in-9 f V Aan/^y' f BRjr ^ ( BR{t ^ e^) \ 

5i?(/. - 67) > 10 (^^J [ g 8 , ,0^3 J 1.1 , 10-^ ) ■ 

(43) 

Therefore, if the rates for the rare tau decays were just below the present experimental 
bounds, the see-saw scenario with universal soft terms at Mx would predict a branching 
ratio for /i — > 67 typically larger than 10~^, in conflict with the present bound BR{fi 
67) < 1.2 X 10"^^. Furthermore, if the observation of both rare tau decays is at the 
reach of present 5-factories, BR{t — > li'-f) have to be larger than ~ 10~^, which would 



■^One loop QED corrections to the electric and magnetic dipole operators reduce the theoretical 
prediction for BR{lj kj) by a factor ^1 — ^ log ^ j [41]. This correction makes the bound 
Eq. (42) a 2-6% stronger for 7713 = 100 - 1000 GeV. 
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imply BR{fi 67) ^ 2 x 10"^'^, barely compatible with the present upper bound from 
MEGA. One can make this argument more quantitative turning Eq. (43) to obtain an 
upper bound on BR{t 67), or alternatively on BR{t fi'-f), from the stringent 
experimental constraint on BR{fi — * 67): 

BR(r e7) < 10^ ( ( E^il^^] ( BR{r^,l) y' 

^ e7j _ iu |^200GeVy V 10 7 V 1-2 x IQ-n J\ 6.8x10-8 J ' 

BR(r ^ U7) < 7 X 10^° (^^) " V ( m^^] ( BRir^e,) \ 
^ /i7j ^ ^ X 1,200 GeVy V 10 J \ 1-2 x 10-^ 7 V 1-1 x 10"^ J ' 

Then, if BR{t — > /i7) ^ 10~®, thus making the observation of r ^ fi'-f accessible to 
present -B-factories, the previous equation would imply that BR{t — > 67) ^ 8 x 10"^, 
which would make the observation of r ^ 67 difficult (analogous conclusions can be 
drawn for BRij /i7))- Therefore, whereas the observation of one rare tau decay at 
present i?-factories is indeed possible, in the see-saw scenario the observation of both 
rare tau decays is unlikely. 

This qualitative discussion shows that in the supersymmetric see-saw model the so 
far negative searches for /i — > 67 have crucial implications for the future searches for 
T ^ fij and r — > 67. Clearly, these implications will become stronger as the MEG 
experiment at PSI improves the bound on /i — s> 67. However, one should bear in 
mind that this is just a qualitative discussion and that the actual impact of the bound 
Eq. (42) on future searches for rare decays depends on the particular point of the SUSY 
parameter space, through the values of and ms". We will see, however, that these 
strong conclusions hold for a wide choice of supersymmetric parameters. 

We have investigated in detail the 'Snowmass Points and Slopes' (SPS) [42], which 
constitute a set of benchmark points in the supersymmetric parameter space aiming 
to describe typical points, but also extreme although well motivated possibilities. We 
have analyzed the six mSUGRA benchmark points, that are defined at Mx = 2 x 10^^ 
GeV by five parameters: the universal scalar mass (mg), gaugino mass (M1/2) and 
trilinear term (Aq), tan/3 and the sign of fi. For each SPS point, the values of these 
parameters are given in Table 2. 

The SPSla and SPSlb points are "typical" mSUGRA points with intermediate and 
relatively high values of tan/5 respectively. Assuming that the neutralino constitutes 
the dominant component of dark matter of the Universe, these two points lie in the 
"bulk" region of the allowed mSUGRA parameter space. On the other hand, the 
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mo (GeV) 


Ml/2 (GeV) 


(GeV) 


tan (5 sign /i 


SPSla 


100 


250 


-100 


10 + 


SPSlb 


200 


400 





30 + 


SPS2 


1450 


300 





10 + 


SPSS 


90 


400 





10 + 


SPS4 


400 


300 





50 + 


SPSS 


150 


300 


-1000 


5 + 



Table 2: Parameters at Mx for the six mSUGRA SPS benchmark points. 

SPS2 point is characterized by heavy squarks and sleptons and fairly light neutralinos, 
charginos and gluinos^. The SPS3 point has a small stau-neutralino mass difference and 
lies in the stau coannihilation region. The SPS4 point has a large tan (3 and lies in the 
"funnel" region. Lastly, the SPS5 point is characterized by a relatively light stop. We 
have also analyzed for completeness the mSUGRA-like scenario SPS6, characterized 
by having non-universal gaugino masses, and defined at the GUT scale by mo = 150 
GeV, M3 = M2 = 300 GeV, Mi = 480 GeV, Aq = 0, tan/3 = 10 and positive fi. 

In Fig. 1 we show the allowed values for BR{t 67) and BRij — > /i7) in the 
super symmetric see-saw model for the mSUGRA benchmark points SPSl-5. The area 
above (to the right of) the dashed line at BR{t — /i7) = 6.8 x 10^*^ {BRij 67) = 
1.1 X 10~^) is excluded by the present experimental bounds on the rare tau decays. On 
the other hand, the area above the diagonal line labeled BR{fi —* e-f) < 1.2 x 10"^^ 
is excluded from the present experimental bound on /i ^ 67, as a consequence of 
Eq. (42). We find remarkable that for all the mSUGRA SPS points the bound Eq. (42) 
excludes values for the branching ratios of the rare tau decays that are otherwise 
allowed by direct searches. Furthermore, if the MEG experiment reaches the sensitivity 
BR{fi — > 67) ~ 10"^^ without observing a positive signal, the region of the parameter 
space excluded by Eq. (42) would enlarge considerably. Finally, in the plots we assumed 
an intermediate decoupling scale, M^aj = 5 x 10^^ GeV. Had we used a larger value 
for Mmaj, the excluded region would also enlarge, as is apparent from Eq. (38)^. 

■*The low energy predictions for this benchmark point are extremely sensitive to the value of the 
top quark mass [43]. In particular, assuming Mt = 175 GeV this point would lie in the "focus point" 
region of the allowed mSUGRA parameter space. However, the most recent measurement of the top 
quark mass at CDF II, Mt = 170.8 ± 2.2(stat.) ± 1.4(syst.) GeV [44], pushes the SPS2 benchmark 
point out of the "focus point" region. 

^Recall that to bring the decay rates to the reach of present and future experiments the neutrino 
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The bound Eq. (42) also has imphcations for future searches for rare tau decays. In 
Fig. 1 we show with a dash-dotted hne the projected sensitivity of present 5-factories 
to rare tau decays {BR{t fi'-f), BR{t 67) ^ 10^^). Therefore, the area shaded in 
green is the region of this parameter space accessible to present 5-factories, that we call 
for definiteness the "observable window of present i?-factories" . Using Eq. (42) we find 
that large regions of the observable window of present 5-factories are excluded from 
the present bound BR{fi —>■ 67) < 1.2 x 10~^^ (light green shaded area). In particular, 
for the benchmark points SPSla, SPSlb and SPSS the region where both t ^ fi'-f 
and r ^ 67 could be discovered at present i?-factories is excluded. Nevertheless, the 
discovery of one of them, either r — ^ /i7 or r ^ 67, is still possible. Let us remind that 
SPSla and SPSlb correspond to "typical" mSUGRA points, and thus this conclusion 
holds for a large region of the parameter space. On the other hand, for SPS2 and SPSS 
the discovery of both rare decays is still possible, although only if their branching ratios 
are close to the experimental sensitivity of present 5-factories^. Finally, for SPS4 the 
present bound BR{fi — 67) < 1.2 x 10~^^ has only little impact for present i?-factories. 
An improvement of the bound on BR{fi 67) by two orders of magnitude, as planned 
by the MEG experiment at PSI, would exclude the possibility of observing both rare 
tau decays at present 5-factories in most mSUGRA parameter space. Hence, would 
present i?-factories refute this expectation, the supersjTumetric see-saw model with 
mSUGRA would have to be abandoned, unless a certain amount of fine tuning is 
accepted. The same conclusions hold for the mSUGRA-like scenario SPS6, which has 
non-universal gaugino masses, as can be realized from Fig. 2. 

We also show for completeness the observable window of the projected supers- 
factories, shown as a yellow shaded area. It is defined by the region between the 
dashed lines indicating the present bounds on the rare tau decays and the dotted 
lines showing the projected bounds {BR{t fij),BR{T fij) > 10"^). The present 
bound on /i ^ 67 practically does not exclude any region of the observable window of 
the projected supers-factories, except for the benchmark point SPSS. On the other 
hand, if the MEG experiment reaches the projected sensitivity BR{fi 67) ~ 10"^'^ 

Yukawa couplings have to be sizable. This typically requires large values of M,„aj in order to produce 
small neutrino masses. 

^One should bear in mind, however, that the SPS2 point is extremely sensitive to the input value 
of the top quark mass. Therefore, the conclusions for this benchmark point should be taken with a 
pinch of salt. 
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without observing a positive signal, again a large portion of the observable window of 
the projected supers-factories would be excluded by Eq. (42). In particular, whereas 
the observation of either r ^ or r ^ 67 will indeed be possible at the projected 
super5-factories, the observation of both would only be possible for the benchmark 
point SPS4 and marginally for SPS2 and SPSS. For the rest of the benchmark points 
analyzed in this paper, SPSla, SPSlb, SPSS and SPS6, the observation of both rare 
decays will not be possible in the supersymmetric see-saw model, unless a certain 
amount of fine tuning is accepted. 

5 Conclusions 

A series of neutrino experiments have demonstrated that all family lepton numbers are 
violated in Nature. Therefore, there is no symmetry reason forbidding the rare lepton 
decays /i — > 67, r — > /i7 or r ^ 67. In this paper we have discussed the implications 
of this observation for the supersymmetric see-saw model. 

We have shown that even in the very special situation where the flavour violation 
vanishes at high-energies in one of the sectors, for instance in the fi—e sector, the flavour 
violation in the other two sectors will induce through two loop radiative corrections a 
small though no n- vanishing flavour violation at low energies in the fi — e sector, that will 
be correlated to the flavour violation in the r— /i and r— e sectors. Barring cancellations, 
this scenario will produce the minimal rate for the rare decay /i — > 67. Therefore, in 
any other scenario the lower bound BR{fi 67) ^ C x BR{t fi'~f)BR{T 67) will 
hold, where C is a constant that depends on supersymmetric parameters. 

We have analyzed the implications of this bound on the possible values of the rates 
of the rare tau decays for the supersymmetric benchmark points SPSl-6. We have found 
that values for BR{t — > fi'-f) and BR{r 67) that are allowed by present experiments 
searching for rare tau decays are forbidden by our bound BR{fi 67) ^ C x BRij — > 
fi'-f) BR{t 67) and the present constraint on BR{fL — > 67) from MEGA. In particular, 
we have found that, for large regions of the Constrained MSSM parameter space, 
present 5-factories could discover either r ^ fi^ or r ^ 67, but not both. We have 
also discussed the implications of the non-observation of the process /i — 67 by the 
MEG experiment at PSI for the future searches for rare tau decays at present B- 
factories and at projected superS-factories. This analysis could also be extended to 
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Figure 1 : Allowed values for the branching ratios of the rare tau decays r — 67 and r — s- /X7 
from present experiments and from the bound BR{fi — > 67) ^ C x BR{t — > fij)BR{T 67) 
for the mSUGRA scenarios SPSl-5. The area in green indicates the observable window of 
present i?-factories, and in yellow the observable window of future super S-factories. Excluded 
regions are shown with light shading, whereas allowed regions are shown with dark shading. 
In the plots it is assumed Mmaj = 5 x 10^'^ GeV. 



SPS6 




Figure 2: As in Fig. 1, but for the mSUGRA-like scenario SPS6. 
more general classes of models. Work along these lines is already in progress [45] . 
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Appendix 

In this appendix we report the full set of one- and two-loop Renormalization Group 
Equations (RGEs) for the parameters of the Minimal Supersymmetric Standard Model 
extended with three right-handed neutrino superfields. Partial results for the Yukawa 
couplings, the gauge couplings and the neutrino mass matrix can be found in [46]. The 
RGEs for the soft SUSY breaking terms and the fi term have been derived particular- 
izing the general formulas in [47] to the supersymmetric see-saw modef. 

The RGE for any supersymmetric or soft-breaking parameter can be schematically 
written as 



d 
It 



-X 



(1) 



(167r2) 



(2) 



(45) 



^All the RGEs are written in the _Di?' scheme. 
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The one- and two-loop (3 functions for the gauge couplings are given by 



(46) 
(47) 



where gi, g2 and are the U{1)y, SU{2)l and SU{'i)c gauge couphngs, respectively, 
and 
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(In these expressions, gi has been normahzed as in S'f/(5).) 

On the other hand, the complete superpotential of the MSSM extended with right- 
handed neutrinos reads, imposing R-parity conservation, 

(49) 

The one- and two-loop (3 functions for these SUSY conserving parameters are: 
P^^l = Y,|Tr(3YrfYi + yM) + 3Y|Yrf + y1y„ - ^gl - ?>gl - ^gl^ , (50) 

(51) 



/5- 



(2) 



Y,| - Tr(9YrfY|Y,Y| + 3Y„Y|Y,Yt + 3Y,Y|YeYt + Y^Y^Y^Yt) 

- YlY„Tr(3Y„Yl + Y.Y^) - 3Y|Y,Tr(3Y,Yi + Y^YJ) 

- 4YiYdYiY, - 2Y t Y„Y t Y„ - 2Y ^ Y^Y^Y, 

IQgl - \g'^TT{Y,Y\) + ^^7?Tr(Y,Yt) + ^^^.^YIy. + kgl + \gl 
5 J 5 5 L 5 



YiY, 



16 



9i 



P^l = Y„|Tr(3Y„Yt + Y.Y^) + SY^Y^ + YIY^ - ^gl - 3gl 



2 13 o 



15 



9i 



(52) 



22 



- Tr(9Y„YlY„Yl + SY^YiY^Y^ + SY^Y^Y^Y^ + Y.YlYeY^) (53) 

- YlYrfTr(3YrfYl + Y^YI) - 3YtY„IV(3Y„Yt + Y.Yt) 

- 4YiY„YiY„ - 2YiY^YiY^ - 2YiY^Y;l Y 

4 



a ^ d ^ u ^ u 



:9l 



X 6 ^ 9 9 

y ^3 + %3^2 



Tr(Y„Yt) + 

136 2 2 



:9l 
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Y^2 + 929l 
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l50 



■^1 



Ye|Tr(3Y,Yj + Y^YJ) + 3Y|Ye + Y^Y, - ?,gl - , 

Y,| - Tr(9Y,Y|Y,Yl + 3Y„Y|Y,Yt + 3YeYiY,Yt + Y^Y^Y^Y^ 
- YtY,Tr(3Y„Yt + Y^Y^) - 3YtY,Tr(3Y,Yi + Y^Yt) 



ayW Y^Y - 2YTY Y^Y - 2YfY Y^Y 

^-'-e-'-e-'-e-'-e ■^-■-i/-'-i/-'-z/-'-// ■^-"-/y-'-iy-'-e-'-i 
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2 2 
-5'\ 



Tr(YrfYi) + ^glTT(Y,Yl) + GglYlY, 



15 A ^22 4 

+ y^2 + 5^2^i + y^i 



(54) 
(55) 



Y, jxr (3Y„Yt + Y.Y^) + 3YtY, + Y+Y^ - 3gl - ^gl | , (56) 

Y, I - Tr (9Y,Y t Y,Y 1 + 3Y,Y1Y,Y t + 3Y,Yt Y.Yt + Y, YIy^YI) (57) 

- Y|YeTr(3YrfYl + Y^YJ) - 3YtY,Tr(3Y„Yt + Y.Y^) 

'^Y j^Y i/Y j^Y 1/ 2YjYeY|Yg 2YeYeYj'^Y^ 

+ \l6gl + ^gllTiiY^Yl) + kgl + ^glWlY, + ^g^YlY, 
L 5 J L 5 J 5 

15 4 9 2 2 207 4 
+ Y92 + -^9.9, + 

= /i|Tr(3Y,Yt + 3Y,Y^ + Y^YI + Y.Yt) - ?,gl - ^-gl^ 

= ^1 - Tr( 9Y„Y,lY„Yt + 9YrfY|Y,Yi + 6Y„YiY,YT 
+ 3YeYeYgYe + 3Y^Yj'^Y^YJ! + 2Y^YgYeYj'! 



(58) 
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- 2Y*YJy:yJ - 2Y*YjY*Yj - 2Y;YjTr(3Y„Y,i + Y.Y^) (61) 
+ ^gfYtYj + 6^22Y*YJ 

- 2Y,YtYeYt - 2Y,YtY,Yt - 2Y,YtTr(3Y„Yt + Y.Yt) 



+ IglY.Yl + GglY^Yt 
5 
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The soft SUSY breaking Lagrangian of the MSSM extended with right-handed 
neutrinos reads: 



soft 



m 



BHuHd — -vpuB-M.ijURj + h.c. 



The /3-functions of the soft gaugino masses are given by 
/3g = 2,,^5(^)M„ , 



(62) 



+ E ^^^(TrfYtAJ-M^TriYtYj) 



x=u,d,e,i/ 



(63) 
(64) 



-b=l 

where Bi^\ B^J and Q-'^-'^.^ were defined in Eq. (48). On the other hand, the (3- 
functions of the soft scalar masses read: 



A 



(1) 



Tr[6(m^^ + m^)YlY„ + 6Ylm2Y„ + 6AtA„ + 2(m^^ + mi)YtY, (65) 
+ 2Ytm2Y, + 2AtA.] - GgllM^l' - ^glW + ^glS , 
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+ 18{ml^ + mj) YtY„YtY„ + 18Ytm^Y„YtY„ 
+ 3(m^^ + mj,^ + m|)YiY„YiY, + SYtm^ Y^Y^Y^ 
+ SYlY^m^YlY, + SYlY^Ylm^Y, 
+ 6(m^^ + mi)YtY,YtY, + GYtm^Y^YtY. 
+ (^k + + mi)YtY,YtYe + Ytm^ Y.YiY, 
+ YtY.miYlYe + YiY.Ylm^Ye + ISA^ A^Y+Y, 
+ ISA^YyY^A^ + 3A|At;YuY„ + 3Y|Yc(A^Au 
+ 3A|YrfYtA« +3Y|ArfAtY„ +6AtA,YtY, +6AtY,YtA, 
-\- A^AgY'jI/Y^^ -\- Y"]! Y'p A;y A;y -I- A^Y'pY^/t A,y Y^pApAjtY^j^ 



e.£-ve-n-j/ ' J/ 



Tr[(m^^ + m^) YtY„ + Yim^ Y„ + A^ A„ 
+ 32^?!^ 2|M3|'Tr[YiYj - M;Tt[YIAu] - M3Tr[AlY„ 



^gl{2\M,\'TT[YlY,,] - A/;Tr[YtA„] - AfiTr[AtYj} + ^g^S' 



+ 33gt\M2\' + fghUm' + |Mi|2 + 3fJ[MiM*]) + ^9tm' 
3 



5 



Tr 



6(m^^ + m^)YiY, + eY^m^Y, + 2(m^^ + mi)YtYe + 2Ytm2Y, 
+ GAlAd + 2A|A, 



6g',m'-lg!m\'-lglS 



2Tr 



+ 18(m^^ + m^)YiY<,YiYrf + ISYjmrfYdYjY, 

+ 3(m^^ + ml^ + m^)YtY„Y|Y, + 3Y,tm2 Y„Y^Y, 

+ 3YtY„m^YiY, + 3YtY„Y|m^Yrf 

+ 6(m^^ + mi)Y|Y,YtYe + eYlm^YeYiY, 



m|)Y^Y,YlYe 



Y,tm?,Y,Y|Ye 



+ YtY.miYtY, + YtY^Ylm^Y, + ISA^ AdY^Y, 
+ ISAlYdYlArf + 3At AnY^Y, + 3YtY,AlA, 
+ SA^Y.YiAd + 3YlA„AiY, + GAlAMYe + 6AtY,Y|Ae 
+ Aj^AjyYjYg + YjI^Y^A^Ag + Ajj^Yj^YjAg + Yt'^A^AgYg 



25 



+ 



32^73' - ^g'l 



Tr[(m|^ + mJ)YiY, + VWaYd + A^Ad] 
+ 32gl[2\Ms\'TT[YlYd] - Af^TrfYU,] - MsTtIMY,]} 

- ^gl[2\M,\'TT[YlY,] - Af*Tr[YiAJ - AfiTr[A|Yj} 

+ ^gf[ Tr[(mi^ + mi)YtY, + Ytm^Y, + AUe] 

+ 2|Mi|2Tr[YtYe] - MiTr[AtYj - Af*Tr[Y|A,]} 

- ^gfS' + 33gtm' + ^glgl{\M,\' + \M,\' + ^M,M;]) 
+ ^g^^\M^\^ + 3gla2 + ^glai , 

(m^ + 2m^jYtY„ + (mj + 2m^jYlY, + 2Ytm2Y„ + 2Yim2Y, (69) 
+ [YtY„ + Y^Y,]m^ + 2AtA„ + 2MA, 

- fgllMsl' - 6glm' - ^^gfm' + ^g^S , 

- (2m^ + 8m|jYtY„YiY„ - 4Ytm2Y,YtY„ - 4YtY„m^YtY„ (70) 

- 4YtY„Ytm2 Y„ - 2YtY„YtY„m| - (2m^ + Sm^^J y1Y<,Y1Y, 

- 4Y|m^Y,Y|Yrf - 4YiY,mjYiY, - 4YjYrfY|m2Y, - 2YiY,YiY,m^ 
[m^ + 4m^J YtY„ + 2Ytm2 Y„ + YtY„m^] Tr(3YtY„ + Y^Y.) 
(m| + 4m^jYjYrf + 2Yim2Y, + Y^Yrfmj] Tr(3YiY, + YIY,) 

- Y,lY„Tr(6mjY,lY„ + 6YTm2 Y„ + 2miYiY, + 2Ytm^Y,) 

- YiY,Tr(6mjYjY^ + QYWaYd + 2miY|Ye + 2Y|m2Ye) 
^•^YyY^AyA^ -|- A^AyY^Y^ -|- Y^A^A^Y^ -|- A^Y^Y^A^ 

- 4{ Y|Y,AiA, + AUdYlY, + YIA^A^Y, + A^Y.YlArf} 

- AtA„Tr[6YtY„ + 2YtY,] - YtY„Tr[6AtA„ + 2AtAJ 

- AtY„Tr[6YtA. + 2YtA,] - YU«Tr[6AtY„ + 2AtY,] 

- AiAdTr[6YjY, + 2YtYe] - YiY,Tr[6AjAd + 2AtAJ 

- AiY,Tr[6YjA, + 2YUe] - YiA,Tr[6AjYrf + 2A|Ye] 

+ ^^^1 (2m^ + 4m^jYtY„ + AYIuiIY^ + 2YiY„mJ + 4AtA, - 4AfiAlY„ 
- 4A^*YtA, + 8\M,\'YlY^ + (mj + 2m^jY|Y, + 2Yjm2Y, 
+ Y^Y.m^ + 2AiA, - 2^/iAiY, - 2Af*YiA, + 4|A^i|2YjY, 
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o 1 98 

+ T^glS' - —gt\M,\^ + W,gl{\M^\' + \M,\^ + a^lM^M*]) 
32 

+ -glgl{\M,\^ + |Mip + 3?[MiM3*]) + ?,?,gl\M^\^ 

9 1 QQ 

+ -glgl{\M2? + iMil^ + 3?[MiM*]) + — <7^|Mip 

16 2 o 2 1 2 



= (2m^ + 4m^jY,Yi + 4Y,mjYl + 2Y,Ylm2 + 4A,A1 (71) 

d 

-^^gl\M,\^ - ^/,\M,\^ + \glS ^ 
(3^^, = - (2m2 + 8m^jY,YiY,Yi - 4Y,m^YiY,Yi - 4Y,Y|m2Y,Yi (72) 

d 

- 4Y,YiY,mjYl - 2YMYdYWd - (2m^ + 4m^^ + 4m^jYrfYtY„Y| 

- 4Y,mjYtY,Yl - 4YrfYtm2 Y.Y^ - 4Y,YtY„mjYi - 2Y,YtY„YK^ 

- (m^ + 4m^jY,Yi + 2Y,mjYl + Y,YW^ Tr(6Y|Yrf + 2YiYe) 

- 4Y,YiTr(3m^Y|Y, + SYjm^Y, + miYlY, + Yim^Ye) 

- 4{ A,A5YrfY| + Y,Y|A,A5 + A^Y^YrfAl + YrfAlA.Y^} 

- 4{ A,AtY„Y^ + Y,YtA.Ai + A.Y^Y^A^ + Y^AtA^Y^} 

- 4A,AiTr(3YiYrf + YJY^) - 4Y,YiTr(3AiA, + A|Ae) 

- 4A,YiTr(3AiY, + A^Y,) - 4Y,AiTr(3YiAd + ^iK) 



+ 



^gl + \gl 



'jnl + 2m^jY<,Yi + 2Y<,mjYi + Y.Ylm^ + 2A,Ai 



+ l2gl\2\M2\''YM - M;A,Y\ - M^Y.M 
+ ^gl{^2\M,\^Y,Yl - M*AaYl - MiY^A^} + ^g^S' 

- ^-fgtm' + ^-^gigum' + + miM;]) + ^-^gtm\' 

16 2 4 2 

(3^^, = {2uil + 4m^jY„Yt + 4Y,m^Yt + 2Y,Ytm2 + 4A„At (73) 

-fgW-%g!m\'-lg!s, 

= - (2m^ + 8m^jY„YtY„Yt - 4Y„mjYiY„Yi - 4Y,Ytm2 Y^Y^ (74) 

- 4Y„YtY„m2^Yt - 2Y„YlY„Y,tm2 - {2inl + Amj,^ + 4m2,jY„YiY,Yt 

- 4Y„m|YlY,Yt - 4Y„Yim2YrfYt - 4Y„YlYrfm|Yt - 2Y„YiYrfYtm2 

27 



>l + 4m^jY„Yi + 2Y„m^Yi + Y^Ylm^] Tr[6YlY„ + 2YtYj 

- 4Y„YtTr[3m^YlY„ + SYim^ Y„ + miYtY, + Ytm^Yj 

4 1^ A„ A,j Yj^ Yjj + Yj^Y^A^A^ + A,jY„Y,jAj^ + YjjAj^A^Yj^ j- 

- 4 { A„AiY,Yt + Y„Y|A,At + A^Y^Y^At + Y.AlA.Yt} 

- 4A„AtTr(3YlY„ + Y^Y.) - 4Y„YlTr(3AtA. + A^A.) 

- 4A„YtTr(3AtY„ + A^Y.) - 4Y„AtTr(3YtA. + Y^A,) 

+ [6^?2 - Idi] {(m^ + 2m^jY,Yl + 2Y„mjY,i + Y^Ytm^ + 2A„A1} 
+ 12gl\2\M2\^YuYl - M^A^YI - M^Y^M 



- -gl^\M^\^YuYl - MlAuYl - MiY^Al} - -glS' 

16 2 16 2 
6 15 

= + (mi + 2m^jYtY, + (mi + 2m^jYtYe + 2Ytm2Y, + 2Ytm2Ye 

+ [YtY, + YtYe]mi + 2At A, + 2At A, 

-6glm'-lgl\M,\'-lglS, 
= - (2mi + 8m2,jYtY,YtY, - 4Ytm2Y,YtY, - 4YtY,miYtY, 

- 4YtY,Ytm2Y, - 2YtY,YtY,mi - (2mi + 8m^jYtY,YtYe 

- 4Ytm2YeYtY, - 4YtY,m? Y^Y, - 4YjY,Ytm!Y, - 2YtY,YjY,n 



e -■- e'-'-'-Ij ^ e ^ e 

[ml + AmljYlY, + 2Ytm2 Y, + Y^Y^mi] Tr(3YtY„ + Y+Y, 

[mi + 4m^jYtYe + 2Ytm2Y, + YlY^mi] Tr(3YiYrf + Y^Y^) 

- YtY,Tr[6mjYiY„ + GY+m^Y^ + 2miYtY, + 2Ytm2Y,] 

- YtYeTr[6mjYiY, + GY^m^Y^ + 2miY|Y, + 2Y|m2Ye] 
4'^Y^Y^A^A^ -|- A^A^Y^Y^ -|- Y^A^A^Y^ -|- A^Y^YJ/Aj^ 

4"|YjYeAjAe ~l~ AjAgYjYe + YlAgAjYg -|- AjYgYjAej' 

- AtA,Tr[6YtY„ + 2YtY,] - YtY,Tr[6AtA. + 2AtAJ 

- AtY,Tr[6YtA« + 2YtA,] - Y^A,Tr[6AiY„ + 2AtY,] 

- AtAeTr[6YlY, + 2YtY,] - YtY,Tr[6AiArf + 2AtA,] 

- AtYeTr[6YlA, + 2Y|Ae] - Y|AeTr[6AiYrf + 2AtYe] 



28 



+ ^gl!^ (mi + 2mljYlYe + 2Ylm2Ye + Y^Yemi + 2AUe 
- 2MiAlYe - 2M;YiAe + 4|Mi|2Y|Ye| 

- ^glS' + 33gt\M,\' + ^glgl{\M,\' + \M,\' + 3?[MiM*]) 

+ ^g^^\M,\' + 3gjcr2 + lgfa, , 
25 5 



(2m2 + 4mi,jYeY| + 4YemiYl + 2YeY|m2 + 4AeA| (77) 

- (2m2 + 8mi,jYeYtYeYl - 4Y,miYiYeYt - 4YeYlm2YeYi (78) 

- 4YeYtYemiY| - 2YeY|YeYtm2 - (2m2 + 4m|^ + 4m^jYeYtY,Y| 

- 4Y,miYtY,Yt - 4Y,Ytm2 Y.YJ - 4Y,YtY,miYl - 2Y,YtY,Ytm^ 

- (m^ + 4mi,jYeYl + 2YemiY| + Y^Ylm^] Tr(6YiY, + 2Y|Ye) 

- 4YeYtTr(3m^YiY, + SY^m^Y, + miY^Y, + Yjm^Y,) 
4-1 AgAlYgYj -|- YgYlAgAl + AgYlYgAl + YeAlAgYj^ 
4"! AgAjiY^Yj + YgYj^Aj/Aj + AgYjjYjyAl + YgAj^Aj^Yl j- 

- 4A,A|Tr(3Y|Y, + Y^Y^) - 4Y,YtTr(3AiA, + AJ Ae) 

- 4AeYtTr(3A|Yrf + A|Y,) - 4YeA|Tr(3YiA, + YUe) 

+ [6gl - Igf] {(m^ + 2mi,jY,Yt + 2Y,miYt + Y^Ylm^ + 2AeA|} 

+ 12gl{^2\M,\'YM - M;A,YI - M^YM} 

12 r "I 

- -gl{2\M^\^YM - M^AeYt - M,YM} 



12 o^, 



2808 4, ,2 12 2 

—-g^^\M,\' + -gfa, 



(2m^ + 4m^jY,Yt + 4Y,miYt + 2Y,Ytm2 + 4 A, A 



t 



(2m2 + 8m^jY,YtY,Yt - 4Y,miYtY,Yt - 4Y,Ytm2Y,Yt 



(79) 
(80) 

- 4Y,YtY,miYt - 2Y,YtY.Ytm^ - (2m^ + 4mi,^ + 4m^jY,YtYeYt 

- 4Y,miY|Y,Yt - 4Y,Ytm2Y,Yt - 4Y,YtYemiYt - 2Y.Y|YeYtm2 



4m^jY,Yt + 2Y,miYt + Y.Y+m^ 



Tr[6YtY„ + 2YtY, 



- 4Y,YtTr[3m2 YlY„ + 3Ytm2Y, + miY^Y, + Y^m^Y, 



29 



4 1^ A;^ AJ/ Y;^ Yj!/ -|- Yj^Y^/A^Aj/ + A^Yj/Yj^Aj^/ + Yj^Ajl/Aj^Yj'/ j- 
4 |Aj^A|YeYji + Yj^YlAgAjji + Aj^YjYgAj'^ + Yj^AjAgY^'i J- 

- 4A,AtTr(3YtY„ + Y^Y.) - 4Y,YtTr(3AtA« + A^A.) 

- 4A,YtTr(3AtY„ + A^Y.) - 4Y,AtlV(3YtA, + Y^ A,) 

+ [6^72' + l9l]{{nil + 2m^jY,Yt + 2Y,miYt + Y.Ytm^ + 2A,At} 
+ 12gl{^2\M2\'Y,Yl - A.f*A,Yt - M^YM} 
+ -gl[2\Mi\^YM - M*A,Yt - AhYMj , 
where we have defined 



S = "^k - + Tr[m^ -ml- 2ml + md + m; 



ej ) 



^1) 



(3m^^ + m^)YtY„ + 4Ytm2 Y„ + (3m^^ - m^)YiY, - 2Yjm2Y, 
(m^^ + mi)YtY, - 2Ytm2Ye + (-m^^ + mi)YtY, 

{m^^ -m^^ -Tr(m2)} + 



3 o 3 o 
2^2 ^Qyi 

16 2^16 , 
y^3 + Y^^i 



;93 



2 , 3 „2 



2^2 3Q^1 



Tr(m^ 



Tr(m^) 



3^3^ + I^^i 



Tr(m^) + -9lTT{nil) , 



a, = ^^7?{3(m|^ + + Tr[mJ + 3mi + Sm^ + 2nil + em^]} 

da = 92 {"^k + "^L + Tr[3mJ + m^]} , 
cr3 = (73^Tr[2m^ + m^ + m^] . 



Finally, the (3 functions for the trilinear and bilinear soft terms are: 

2 



A,|Tr(3Y,Y^ + Y,Yt) + SY^Y, + YtY„ - 
Y JTr(6ArfY| + 2AeYi) + 4YiA, + 2YUn 



16 2 2 2 

y^3 -3^2 - Y^^i K 



^2) 



32 



14 



—9iM3 + 6g'^M2 + Y^9{Mi 



(2) 



A,| - Tr(9Y,Y|Y,Y| + 3Y„Y|Y,Yt + 3Y,YtYeYt + Y^YlY^Yt) 

- YtY„Tr(3Y„Yt + Y^Y^) - 5Y|Y,Tr(3Y,Yi + Y^YI) 

- 6YiY,YiY, - 2Yt Y„YtY„ - 4YtY„YiY, 



(83) 



+ 



16^71 - ^gl\TTiY,Yl) + ^gfTriYM) + ^91yIY^ + 



12gl + -g-^ 



6 
5' 



YiY. 



16 



,2„2 



15 



-"4io22i"22i '-^4i 22 
- Y^3 + 8^3^2 + q939i + + 929l 



287 
"90 



■9i 



30 



Y,<| - 2Tr( 18A,Y|YrfY| + 3A„Y|Y,Yi + 3A,YiY„Y| 

AeYjI^YjyYl) 



6 AeY^YgYg 



4 



A YTY YT 



- YtY„Tr(6A„Yt + 2A,Yt) - 6YlY,Tr(3A,Y| + A^YI) 

- YtA„Tr(6Y„Yt + 2Y,Yt) - 4YlA,Tr(3Y,Y| + Y.YI) 

- 6YiY,Y| A, - 8Yi ArfY^Y, - 4YiA„YlY„ 



4YiY„Y2 A„ - 4Y2 A.Y^Y, - 2YiY„Y^ A, 



+ 
+ 



4 2 



6g. 



19 8 
Tr(A,Yl) + -glTiiAM) + -g^YlA^ 



4 1 12 
32glM, - -(7?MijTr(Y,Yi) - -glAhTT{YM) 



l2glM2 + -glM^ Y\Y, - -g^M^YlY^ 

5 J 5 

64 16 
+ -giMs - IQgjgliM, + M^) - -glgl{M, + M{] 

574 

- 30(72'M2 - 2glgl{M2 + Mi) - -^gtKh 

45 

= A„|Tr(3Y„Yt + Y.Yt) + 5YtY„ + Y^Y, - ^^gl - 'igl - ^^(^i }- 

Y„|Tr(6A„Yl + 2A,Yt) + 4YtA« + 2YIK^ 

32 26 
+ —glM^ + ^glM2 + —g\M^ 
3 15 

/3g = A„| - Tr(9Y,YtY„Yt + 3Y„YlY,Yt + 3Y,YtY,Yt + Y.^t 

- YiY,Tr(3Y,Y| + Y^YI) - 5YtY,Tr(3Y„Yt + Y.Y^ 

- 6YlY„Yt Y„ - 2YlY,YiY, - ^YIY^YIY^ 

+ \\^gl + ^(7i'lTr(Y,Yl) + 12^?22YiY, + ^^^^Y^Y,, 
L 5 J 5 



Y Y^Y Y^\ 



16 A 9 9 136 9 9 

+ %3^2 + -^9z9i 



15 



^2 + aWi 



2743 



Y„| - 2Tr( 18A„YtY„Yt + 3A„Y^Y,Yt + 3A,YtY„Y, 

+ A,Y 



- 6YtY„Tr(3A„Yt + A.Y^) - Y^Y,Tr(6A,Yi + 2A,Yt) 

- 4Yt A«Tr(3Y„Yt + Y.Y^) - Y5AdTr(6Y,Yi + 2YeY|) 

- 6YtY„Yt A„ - 8Yt A^YlY, - 4Y|Y,YiA, 

- AY\KdY\Yd - 2YiY,Yt A„ - 4YiAdYtY„ 



^5) 



31 



Tr(A„Yt) 



5 



4 



u u 



32glM3 + -(??MiJ Tr(Y„Yt) - \ i2glM2 + -g^M, 

- ^gfM.YlY, + ^g^M, - 16glgl{M, + M^) - '^glgl{M^ + Mi 

- 3052'^2 - 2glgl{M2 + Mi) - ^^7iMi 



Ae|Tr(3Y,Yi + Y^YJ) + SY+Y^ + Y^Y, - ?>gl - ^-gl^ + 
Y,|Tr(6A,Yl + 2A,Y|) + 4Y|A, + 2YtA, + QglM^ + y^i^ij , 
Aej - Tr(9Y,YiY,Yi + 3Y,YiY,Yl + 3YeYlYeY| + Y,Y|YeY^) 



YtY,Tr(3Y„Yt + Y.Y^) - 5YtY,Tr(3Y,Yi + Y,Y, 
6"V"|'Vg'Vg"Ve 2'Vj!/'Vj/"V/|/"V// 4'Vii/'Viy'Vg"Ve 



+ 



i%3 - ^gi 



Tr(Y,Yi) + lgli:i{Y,Y, 
5 



12^2^ - Igl 



YlY, 



- 2Tr( 18A,YiY,Yl + 3A,YlYrfYt + 3ArfYiY„Yi 

+ GAgYjYeYj + A^YjYgY^ + A.^YyY yY]^ 

- YtY,Tr(6A„Yt + 2A,Yt) - 6YtY,Tr(3ArfYj + A^YI) 

- Yt A,Tr(6Y„Yt + 2Y,Yt) - 4YiA,Tr(3Y,Yi + Y^YI) 

- GYlYeYl A, - SYl AeYlYe - 4Yt A^Y^Y, 



+ 



V ^ V ^ v-^v 

- pi 



V ^ V ^ e-^e 



12 



Tr(A,Yi) + —gfliiAM) 



Qgl + Igl 



(86) 



(87) 



4 1 12 

32glM, - -glM, Tt{Y,YI) - -(^?AfiTr(Y,Yt) - Ug^M^YlY, 

5 J 5 

SOg^M^ - ^glgliM, + M,) - 54gtM,\ , 



A,|Tr(3Y,Yt + Y.Y^) + SY^Y, + YJY, - 3gl - ^g{^ + 
Y,|Tr(6A„Yi + 2A,Y+) + 4YtA, + 2YUe + 6glM2 + 



32 



Tr(9Y„YlY„Yi + SY^Y^Y^yI + SY^Y^Y^Y] + Y^Y^Y^Yi) 
YtYeTr(3YrfY| + Y^Yl) - 5YtY,Tr(3Y„Yt + Y.Yt) 
fiY^Y Y^Y - 2Y^Y Y"'"Y - 4Y"'"Y Y^Y 

^■'-v-'-i'-'-i'-'-i' '^■^e-^e^e^e ^^e^e-^u-^u 



+ 



-'-^4,922 
Y^2 + -^9291 



Tr(Y„Yt 
207 



12 



-9i 



2Tr( 18A„YtY„Yt + SA^YjYrfYi + SArfY^Y^Yi 



6 A^Y^Y^Y^ 



tv vt _i_ A v'^Y Y"'" -I- A Y'^Y Y'^l 



- 6YtY,Tr(3A„Yi + A^Y^) - Y|YeTr(6A,Yi + 2AeY| 



4YtA,Tr(3Y„Yi 



Y,Y^) 



YlAeTr(6Y,Yi + 2Y,Yt 



GYtY^YiA, - 8Yt A.YiY, - 4Y|Y,YlAe - 4Y|AeYtYe 
2Y|YeYtA, - 4Y|AeYtY, + [?>2gl + ^(7?] Tr(A„Yt) 



+ 



6^72' + \9l 

12 

l2glM2 + —glK'h 
5 



^9) 



Yt A, + ^glYlA, - [?>2glM, + \9lMi] Tr(Y„Yt) 
1 2 

YlY, - -glM.YlYe 



30g',Ah - YdhfiM, + Ml) - ^9tMi 



1 r ' 



i?|Tr(3Y„Yt + SYrfYi + Y.Yt + Y^Yt) - 3gl - } + 
/i|Tr(6A„Yl + 6A,Yi + 2A,Yt + 2AeY|) + 6glM2 + 
^1 - Tr( 9Y„YtY,Yt + 9YdYlYdYl + 6Y„YlY,Yt 



+ 



l6^73' + ^9! 



3YMyM + SYeYlYeYl + 2Y,YIY,YT) 

t 



15 4 9 2 2 207 4 



Tr(Y„Y2) + 

9t 



16^71 - -9! 



6 



(90) 



(91) 



Tr(YrfYl) + -^7?IV(YeYt) 



/i| - 4Tr( 9A„YtY„Yt + 9A,YiY,Yi + SA.YlY.Y.t + SA^Y^Y^Yi 
+ 3 Ajy Yjy Yjy Yjy + 3AeY|YeYl + Aj^Y^YeYT + AeYTY,,Y|) 



+ 



32(71 + -gf 



Tr(A,Yi) 



32(73 - ^9! 



12 

Tr(ArfYl) + y (/?Tr(AeY|) 



32(732M3 + -(7iMiJ Tr(Y„Yl) - [32(?2M3 - -glM,\ Tr(Y,Yi) 
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12 



414 



(1) 

Bm 



(2) 



Bm 
Bm 



+ M 



4Y:aJ 



+ 



2Y,Yt 



B 



M 



4A,Yi 



M, 



2Y*Yf Y:yJ - 2Y;yJy:yJ - 2Y:YjTr(3Y„Yt + Y.Y^) 



(92) 
(93) 



5 



+ 



M 



- 4Y:ArY:Yj - 4Y;aJy;yJ - 4Y;YjTr(3A„Yl + A.Y^) 

- 4Y*YrY: Aj - 4Y;yJy; Aj - 4Y; AjTr(3Y„Yt + Y.Y^) 

12 12 
+ —glYlAl + l2glYtAl - —gfA^YtY^ - Ug^M^YtYj 

- 2Y,YtY,Yt - 2Y,YtY,Yt - 2Y,YtTr(3Y,Yt + Y.Yt) 
6 



+ 



:glY,Yl + QglY^Yl 



B 



- 4Y,YUeY^ - 4Y,YiA,Yi - 4Y,Y^Tr(3A„Yl + A.Y^) 

- 4A,YtY,Yt - 4A,YtY,Yt - 4A,YtTr(3Y„Yt + Y.Yt) 

+ ^-^glA^Yl + l2glA,Yl - ^^glM^Y^ - Ug^M^Y^ 

5 5 



M . 
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